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o 

We discuss the electromagnetic properties of both a charged free particle, and a charged 
particle bounded by an harmonic potential, within collapse models. By choosing a particu- 
, larly simple — yet physically relevant — collapse model, and under only the dipole approxima- 

tion, we are able to solve the equation of motion exactly. In this way, both the finite time 
and large time behavior can be analyzed accurately. We discovered new features, which did 
Oh' not appear in previous works on the same subject. Since, so far, the spontaneous photon 

emission process places the strongest upper bounds on the collapse parameters, our results 
call for a further analysis of this process for those atomic systems which can be employed in 



experimental tests of collapse models, as well as of quantum mechanics. 
PACS numbers: 03.65.Ta, 42.50.Lc, 42.50.Ct 
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Models of spontaneous wave function collapse provide a description of quantum 

(and classical) phenomena, which is free of the much debated measurement problem affecting 
$_i ■ the standard quantum theory. This is achieved by modifying the Schrodinger equation, adding 



nonlinear and stochastic terms which reproduce, at a suitable scale, the process of wave function 
collapse. 

By modifying the Schrodinger equation, these models make predictions which differ from those 
of standard quantum mechanics. It is interesting and important to analyze such differences, not 
only for a better understanding of these models, but also for deciding which experimental setups 
are more convenient, in order to test them against quantum mechanics. Such experiments, needless 
to say, would represent important tests also for the quantum theory itself 
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It has emerged from the work reported in 



10 



111 ], that the electromagnetic properties 



of matter place, so far, the strongest upper bound on the collapse frequency A GRW of the GRW 
(Ghirardi-Rimini- Weber) model [l[l (or, equivalently, the parameter 7 of the CSL 1 (Continuous 
Spontaneous Localization) model 3J). More specifically, it has been proven that charged particles 
spontaneously emit radiation, as a consequence of the interaction with the collapsing field, also 
when according to standard quantum mechanics no radiation should be emitted; the radiation 



spectrum has been computed both for a free charged particle [10] and for an hydrogenic atom [111 ]. 
The theoretical spectrum has been compared with available experimental data, placing an upper 
bound 111] of only 6 orders of magnitude away from the standard CSL value 2 7 = 10~ 30 cm 3 s . 



Note that more direct experiments of the superposition principle of quantum mechanics, such as 

nr 



diffraction experiments with macro-molecules [121. Il3j]. place a much weaker upper bound, which is 
13 orders of magnitude away from the standard CSL value 141 ] . These figures show that analyzing 
the electromagnetic properties of matter within collapse models is particularly relevant, not only 
per se but also in view of future experimental tests. 

The above mentioned analysis has been carried out to first order in perturbation theory, using 
the CSL model. Goal of this work is to deepen our understanding of the process of spontaneous 
photon emission from charged particles. We will do it by using, in place of the CSL model, the 

and 

we will work under the dipole approximation. These assumptions will allows to solve the equations 
of motion exactly: we will derive an exact formula for the spectrum of the emitted radiation, valid 



3, 



111 ]. As we will see, new 



to all orders, and we will compare it with the formulas obtained in 
features will emerge, previously not discussed. 

The QMUPL model of spontaneous wave function collapse applies to systems of distinguishable 
non-relativistic particles. The one particle equation, which is sufficient for the purposes of this 



1 The GRW and CSL models are the two most popular models of spontaneous wave function collapse. Their main 
difference is that the first assumes the collapses to occur at discrete random times, through a jump process, while 
the second assumes the collapse to occur continuously, through a diffusion process. 

2 This value for the collapse parameter 7 has been chosen in such a way that for a single constituent — in which case 
the GRW and CSL model coincide at the statistical level — the reduction occurs with the rate Agrw — 2.2 x f0 _17 s _1 
of the GRW model The relation between the two constants is: Agrw = jia/Aiv) 3 ^ 2 , with a = 10 10 cm -2 [jj. 
This choice implies that the two models behave similarly, although important differences arise, due to the fact 
that the GRW model refers only to systems of distinguishable particles, while the CLS model takes into account 
also identical particles. The numerical value for Agrw has been originally chosen in such a way to ensure that 
superpositions of macroscopic objects (containing roughly an Avogadro's number of constituents) localized within 
the perception time of a human being, while microscopic systems retain all their quantum properties 
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paper, reads: 



dipt 



-~Hdt + VA(q - (q) t ) dW t - ^(q - <q>t) 2 dt 



(1) 



its generalization to a many particle system being straightforward. In the above equation, H is 
the standard quantum Hamiltonian of the particle, q its position operator, Wj three independent 
standard Wiener processes defined on a probability space (f2,JF, P), and A is a positive constant 
which sets the strength of the collapse mechanics. The physical content of the above equation is 
the following. The first term on the right-hand-side gives the usual unitary evolution, driven by 
the Hamiltonian H. The second and third terms cause the collapse of the wave function towards 
a state which is localized in space, being it driven by the position operator. More specifically (but 
not entirely correctly; for a discussion see [13] ) , the third term localizes the wave function — being it 
the more negative, the greater the difference |q— (q)t| — while the second term, which contains the 
random process W(, ensures that the collapse occurs randomly and in agreement with the Born 
probability rule. The structure (q— (q)t) ensures that the wave function remains normalized, even 
if the dynamics is not unitary anymore. 

This model, in spite of its simplicity, is particularly relevant because, in an appropriate limit 
which we shall now briefly discuss, it reduces at the statistical level to the more familiar GRW 
model (thus also to the CSL model, as long as the particles are distinguishable). Here again we 
limit our consideration to one single particle. The master equation describing the time evolution 
of the statistical ensemble pt = ^\\ipt) (ipt\], where ipt evolves according to Ec^Jl]), has a Lindblad 
form, which in the position representation (where pt(x,y) = (x.\pt\y)) reads [16J: 

^Pt(x,y) = ~ [H,p t (x,y)] - ^(x-y) 2 /9t(x,y). (2) 



On the other hand, the one-particle GRW (and CSL) master equation reads 



j t Pt(x,y) = -^[^,p t (x, y )]-A GRW [l-e-f( x -y) 2 ] Pt (x, y ); (3) 

for the relation between the constant A GRW characterizing the GRW model and the constant 7 
defining the strength of the collapse process in the CSL model, see the previous footnote. The 
second parameter (a) defines a correlation length rc = \j\fa ~ 1CP 5 cm, above which spatial 
superpositions are reduced. 

Let us now consider situations where, for all values of x and y such that the density matrix 
pt(x,y) is appreciably different from 0, one has: |x — y| <C rc- We call this the small distances 
assumption. This is the case if the physical system is localized well below rc, as it happens e.g. for 
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atoms in matter. In this case, it makes sense to take, in Eq. ([3]), the limit a — > and A GRW — ► do, 
while keeping the product A GHW a constant. Then, Eq. ([3]) reduces to ([!]), with the identification: 

> "A GRW a 3 / 2 7 

A " ~ 16^" (4) 

Accordingly, the QMUPL model represents, at the statistical level, a good approximation of the 
GRW models, for those systems which are well localized with respect to the correlation length ro- 



ll. MOTION OF A CHARGED PARTICLE INTERACTING WITH THE 
ELECTROMAGNETIC FIELD, BOUNDED BY A LINEAR FORCE, SUBJECT TO 

COLLAPSE IN SPACE 

In this section we explicitly solve, under only the dipole approximation, the equations of motion 
for a non-relativistic charged particle interacting with the second quantized electromagnetic field. 
The particle is bounded by an harmonic potential — the limit case of a free particle will be also 
discussed — and is subject to spontaneous collapses in space according to the QMUPL model. 

Eq. (pQ) is nonlinear, but it can be appropriately reduced to a linear (though not norm-preserving) 
equation trough a standard procedure [lg]. Of course, non-linearity is not canceled; it appears back 
again when the statistical properties (through a change of measure) are computed. However, since 
we are ultimately concerned only with physical quantities of the type Ep[('0t|O|V ; t)], where O is 
any suitable self-adjoint operator, we can use the following mathematical property. Consider the 
class of SDEs: 



d4 



l -Hdt + ^A(Cq - CR<q>t W - ^A(|C| V - 2CC R q(q)t + £(<l)t)dt 



4, (5) 



where £ is a complex phase, and £ R its real part. A straightforward application of Ito calculus 
allows to prove that Ep[(^|0|^)] is independent of £, in spite of the fact that Eq. ([5]) describes 
completely different evolutions for the wave function, for different values of £. In particular, when 
C = 1, Eq. ([5]) coincides with the QMUPL collapse equation (pQ). On the other hand, when £ = i, 
Eq. © reduces to the simpler SDE: 



dipt 



jHdt + iVXqdW t - ^Aq 2 cft 



ipt, (6) 



where all the non-linear terms have disappeared. Such an equation of course does not lead to the 
collapse of the wave function, since it describes a linear and unitary 3 , though stochastic, evolution. 



3 The third term of Eq. ([6]) is an "Ito term", which disappears from the solution of the equation. For this reason 
the evolution is unitary, even if apparently it does not look so. 
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Nevertheless it is as good as Eq. ([I]) for computing averages quantities. The advantage is that its 
linearity and unitarity make calculations easier. 

Eq. ([6]) has to be understood in the Ito sense. We will solve the corresponding Stratonovich 
equation, where the stochastic differential dW t can be interpreted as the increment of a white noise 
w(i): 

ihj^t = [H- A/Aft qw(t)] Vt- (7) 

This is a standard Schrodinger equation with a random potential depending on the position q 
of the particle. Note that the last term of Eq. ([6]) has disappeared in going from the Ito to the 
Stratonovich formulation of the SDE. Actually, one can be more general and assume that w(i) 
represent three Gaussian noises with zero mean and a general correlation function, without having 
to change the mathematical formalism. However, this goes beyond the scope of the present analysis, 
so we will keep assuming that w(t) are white noises. 

Coming back to our physical system, the standard Hamiltonian H is: 

H = ^-(p- e A) 2 + i K q 2 + i e o/^[E 2 + c 2 B 2 ], (8) 

where mo is the bare mass of the particle, k is the force constant of the harmonic term, A is the 
vector potential, E and B the electric and magnetic fields respectively, e is the electric charge, c 
the speed of light and eo the vacuum permittivity. Throughout this section, we use the gauge: 
V • A = and V = 0, where V is the electromagnetic scalar potential 4 . 
The plane wave decomposition of the vector potential A reads: 

A (x) = My I -$L= 6 k „ L, e ik ' x + 4, e-H , (9) 

where uJk = ck (k = |k|) is the frequency corresponding to the wave vector k; (p = 1,2) are 
the linear polarization vectors and aj. , a k/J are the creation and annihilation operators satisfying 
the standard commutation relations: 

KXvl = v 5(3) ( k - k ')- ( 10 ) 



4 Note that we are assuming that the spontaneous collapse process occurs only for the particle, not for the electro- 
magnetic field. The reason is that, so far, collapse models have been considered only for massive particles, being the 
localization of their wave function sufficient for solving the measurement problem of quantum mechanics. However, 
in a more speculative scenario, e.g. where the collapse mechanism is linked to gravitational phenomena, one could 
assume that also the photons' wave function undergoes a spontaneous localization process. 
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Up to now the model is exact, but not exactly solvable. To further proceed in the analysis, we 
make the dipole approximation e lkx ~ 1, which holds as long as the wave-length of the electro- 
magnetic radiation is much larger than the typical size of an atom. Note that this assumption is 
compatible with the small distances assumption discussed in the previous section. The resulting 
model turns out to be ultraviolet divergent: we cure this problem by introducing a form factor 
<?(k), corresponding to the Fourier transform of the charge distribution (normalized to unity): 

d 3 k 



9W 



p{v)e 



- tkr 



Under these approximations, the vector potential (J9j) becomes: 



d 6 r p(r) 



1. 



(11) 



A(x) 



'2u k 



(12) 



and the factor (2ir) 3//2 as been included in the definition of g(k). In this way, in the point-particle 
limit, g(k) -> l/y/pOj*. 

Since the total Hamiltonian in Eq. ([7]) is a standard — though stochastic — Hamiltonian, one can 
conveniently work in the Heisenberg picture. The equations of motions for the position q(i) of the 
particle and of the conjugate momentum p(t) can be immediately derived: 



-j?- = -sq + \/Xhw(t), 



(13) 



dq _P e_ A 

dt 



m 



m 



while the equation of motion for the electromagnetic-field operator (t) is: 



dt 



le 



9(k) 



•E 



d 6 k 



(14) 



a k'fj,' + a k'/i' 



(15) 



the equation for (t) can be obtained from the previous one by taking the hermitian conjugate. 
The above set of coupled linear differential equations can be conveniently solved with the help of 
the Laplace transform; the equations for the transformed variables (which are denoted by a tilde) 
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read: 



zp(z) - p(0) = — nq(z) + hVXw(z), 



(16) 



m m V e 



zq(*)-q(0) = /^-y^5(k)6 kM [4(z)+a kM (z)], (17) 



« 5 L(*) - a L»(°) = 



t 



1 



k^ 



g{k) e k)1 ■ [zq{z) - q(0)] , 



(18) 



z %» - a k M ^°) = ~ iuJk %» + 7y^= 5(k) e kM • [«q(«) - q(0)] , 



(19) 



where z is the transformed time. The above set now represents a system of coupled algebraic 
equations, which can be solved in a standard way. The calculation is long but straightforward; 
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transforming back to the original variables, one obtains: 



q(t) = [l-«Fi(t)]q(0) + F (t)p(0) 



e 



/' 



/2uj k 



G+CMK^OHG^M)^) 



Jo 



(20) 



p(t) = - K [t-KF 2 (t)]q(0) + [l-«Fi(t)]p(0) 



+«e\/ — > 



d 6 k 



~k/x 



G+(MK M (0)+Go(M)<(0) 



/ ds [1 - as Fi (t - s)] w(s) 

JO 



(21) 



e^*4 M (0) - ^^e kM .[Gr(M)p(0)-KG o -(M)q 



| *e 2 g(k) 



e v 7 ^ kAt ^ 



V 2uJ k> 



G- (k, k', t) a k , , (0) + GI (k, k', t) al , (0) 



-le 



(22) 



e-^*a k/i (0) + ^J^e^ - [Gf(k,t)p(0) - KG+(k,t) q (0)} 



ie 2 g(k) 



eo \Z2uJk 



G+ (fe, k', t) a k>/ (0) + Gt (k, k', t) a{^, (0) 



-ie 



(23) 



In the previous formulas, we have introduced the following functions: 



F n (t) = J 

Gt(k,t) 
G%{k,k't) 



dz 

r 27ri z n H(z) 
dz 



z n e zt 



r 27ri (z ± iu}f.)H(z) ' 



n = 0,1,2, 
11 = 0,1, 



z 2 e rf 



(24) 
(25) 
(26) 



/r 27ri (z ± iuk)(z ± iuJk')H(z) ' 

in the third expression, the upper ± refers to the first parenthesis, while the lower one refers to 
the second parenthesis. In all the above formulas, according to the theory of Laplace transform, 
the contour T must be a line parallel to the imaginary axis, lying to the right of all singularities of 
the integrand. The above solutions should be compared with those obtained in [l8,[l^], where the 



collapse process was not taken into account: a part from a marginal calculational mistake in [1 
in the evolution of p(i), the two results agree when A is set to zero in Eqs. (|20j) - (|23l) . 
The function H{z) is defined as follows: 



H(z) = k + z 2 



87re2 r J; ma fc2 
rao + - — / dk g(k) 



(27) 



3e 7 "" JV " y ^ 2 + w 2 _ 
(from now on we assume the form factor to depend only on the modulus k of k). This is a 
crucial quantity, as through formulas (|24p - (|26p it determines the time evolution of all physical 
quantities. It depends on the form factor g(k): simply removing it, would make the integral 
ultraviolet divergent. To overcome the problem, we apply a renormalization procedure. The 
quantity within square brackets in (|27p can be rewritten as follows: 



87re 2 f°° 2 k 2 ( 4 \ 8^e 2 2 f°° g{k) 2 



m + — — / dkg{k) — ^ = [ m Q + -m r ) - »z / dk -7— — (28) 

3e Jo z + u k \ 3 / 3e o c Jo 2 + u£ 

where m r is the electrostatic mass: 



m r := ^ [ dW P -^l ^ ^ I d'kg(k) 2 . (29) 



87reoc^ J |r — r'l e^c 2 



When ^(k) — > \j \J (27r) 3 , m r diverges. We apply the classical renorma' 



20 



izati on p rocedure 5 for a 
221 ] (which is valid 



21 



non-relativistic charged particle coupled to the electromagnetic field 
both as a classical calculation and as a Heisenberg picture, quantum mechanical one, like in our 
case). According to it, as m r — > +oo in the point-particle limit, one assumes that mo — * — oo, in 
such a way that m := rriQ + (4/3)m r remains finite. This is assumed to be the renormalized mass. 

The last term of (|28p instead remains finite in the limit, the integral can be evaluated, and 
H(z) takes the well-behaved expression: 

e 2 

H(z) = K + z 2 [m-Pz], (3 = ;r ~ 5.71 x 1(T 54 Kg s. (30) 

Note that (3 is precisely the coefficient in front of the Abraham-Lorentz force, which is responsible 
for the runaway behavior of the corresponding Abraham-Lorentz equation, as we shall soon see. 
H(z) is a polynomial of third degree, whose zeros can be found by the standard Cardan method. 



One solution is real and two are complex conjugate. Let ujq := \l n/m be the frequency of the 
oscillator. By assuming uj <. 2m/v / 27/3 ^ 6.14 x 10 22 s^ 1 for an electron (huj <C 4.04 x 10 4 KeV), 
their approximate value is (see Appendix A): 



One can notice that the collapse terms do not enter the following equations, thus the renormalization procedure 
applies like in standard cases. 
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Given the above results, the functions F n (t) and G^(k,t) defined in (f24"j) and ([25]) . which are the 
only ones we will use in the subsequent analysis, become: 



Fn(t) = 



zee 



Zft 



H(z) 



+ < 



z=z e 



G±(M) = £ 



z n e Zet 



-* (ze ± ik) 



z- z t 
H(z) 



+ 



n 


= 0, 


uJq 2 n 


= 1, 


tujQ 2 n 


= 2, 


(±ik) n 





H(±ik) 



(32) 



(33) 



The term in (|32p and (|33p with t = 1 diverges exponentially, since z\ is positive. As we have 
anticipated, this is a manifestation of the runaway behavior of the Abraham-Lorentz equation 



21 



22 



20, 



23J. In particular, in the free particle case (ujq = 0), the coefficient z\ corresponds to the 



rate of exponential growth of the acceleration, as discussed in textbooks. This problem is still open, 
and we pragmatically dismiss it by ignoring, in the subsequent formulas, all terms with I = 1. 



A. The spectrum of the spontaneously emitted radiation 

We are now in a position to compute the spectrum of the radiation spontaneously emitted by 
the particle, due to the interaction with the noise. Let iV k (t) := a k „(i) be the density 
of photons of wave vector k and polarization [i. Let \4>) := |^ ho ) |Q) be the initial state of the 
system, where IVO is the initial state of the harmonic oscillator and \Q) is the vacuum state for 
the electromagnetic field. Let finally S*(k, fi, t) := Ep[(0|A 7 k/i (i)|(^)], be the spectrum of the emitted 
radiation, averaged over the noise. By inspecting Eqs. (|22l) and ()23|) . one can notice that all terms 
of -/V k (i) containing either a kA1 (0) or a kfl (0) give a zero contribution, when averaged with respect 
to the vacuum state, while all terms containing w(i) give a zero contribution, when the stochastic 
average is taken. Accordingly, 5(k, /x, t) is the sum of two terms: 

S(k,i2,t) = 5 qm (k,/x,t) + S ml (k,(j,,t), (34) 

where 5 qm (k, fj,,t) is the standard quantum formula, while S col (k, fi,t) is the contribution due to 
the noise. We are interested in computing this second term, which reads: 

\fj p 2 i ft 

S col (k,(x,t) = — 3 / dsG^(k,t-s)Gt(k,t-s), (35) 

167r d eo uJk Jo 

This is the main formula. In the next section, we will apply it to the two interesting cases of a free 
particle (u>q = 0) and of a bounded particle (wq 0). 
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III. THE FREE PARTICLE 



The free particle evolution can be deduced from the previous formulas by taking the limit 
cjo — * 0. However it turns out to be easier to re-do the calculation, starting from Eqs. (113j) - (115p 
with k = 0. The final result is: 



q(t) = q(0) + Fo(t)p(0)- e yj]r J d 3 k^Le kfl 

+\ r \h [ dsF (t - s)w(s), 
Jo 

p(t) = p(0) + y/\h [ dsw(s), 
Jo 



G?(M) a^(0) + GnM)<(0) 



(36) 



(37) 



, 5(k' 



v r 3 



(£T_(fc, fc', i) a kV , (0) + GI (A;, fc', i) a kV (0) 




eo V^k kAt 



(is G l (k, t — s)w(s), 



(38) 



with: 



"^« ku (0) + i4=G-+(M)e k/1 .p(0) 



ie 2 p(k) 



eo 



2w fc k ^ ^J kM 



Gjt(fc, fc', t) a k> , (0) + Gt (k, k', t) a k , u , (0) 



. fhX g(k) f* , - , , , . 

~* e V 7^= e ku ' / dsG{(k,t- s)w(s), 



*b(*) 
G±(fc,t) 



A + HL e rr*/0 



m 



re 



mt//3 



+ 



muo K uo K (m±iuj K ) (m / (3)[(m / (3) ± iu K ] ' 



,^iw it 



+ 



e mt/f3 



+ 1 



(39) 



(40) 
(41) 

(42) 



[(m/P)±iw K ][(m/p)±iu K >] 

In the last expression, the upper ± refers to the sign in front of each u K , while the lower db refers 
to the sign in front of each u K >. Once again, in all above formulas we have a run-away behavior, as 
consequence of the renormalization procedure. In the subsequent analysis, we neglect such terms. 
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There are two quantities which are of particular interest, in order to understand the behavior 
of the free charged particle under the influence of the collapsing field: the evolution of the mean 
kinetic energy, and the spectrum of the emitted radiation. We shall now discuss both of them. 

A. The mean free kinetic energy 

The mean kinetic energy of the particle is given by: 

E meaa (t) = imE[(0|q(t) 2 |0)]. (43) 

From Eqs. (|36|) and (|40p we have: 

+^ r«fow(a). (44) 



m JO 

By taking as initial state \4>) = \tp ficc )\£l) , as in the previous section, and after differentiating over 
time, one obtains the following expression: 

d _ 3 Xh 2 _ 3 A ghw q grw ?i 2 . . 

dt mean ~ ~ 4 m ' 1 j 

which corresponds to the standard GRW formula [l|. We have then a very interesting result: 
in spite the fact that — as we shall see in the next subsection — the particle emits radiation at a 
constant rate, its mean kinetic energy increases steadily in time as if the particle were neutral. In 
other words, the noise drives enough energy into the particle both to increase its kinetic energy 
and to make it radiate. This is a consequence of the fact that the collapse terms contain only 
the position operator q, due to which w< acts like an infinite temperature noise; this feature has 
been first pointed out in [24|. In the same reference, it was shown that a term proportional to the 
momentum operator acts like a dissipative term, thanks to which the mean energy thermalizes to a 
finite value, associated to a temperature which can be considered as the tem per ature of the noise. 
This is similar to what happens in the theory of quantum Brownian motion 25. I^. I^. and more 
generally in the theory open quantum systems, which does not come as a surprise, since collapse 
models and open quantum systems rely on similar master equations. 

The above results can be read in two different way. On a more conservative level, one can 
accept this steady energy increase as a feature of the model, as long as it does not violate known 
experimental data. On a more speculative level, it suggests that the coupling between the noise 
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and the wave function should be modified in order for the total energy (energy of the noise, plus 
kinetic energy of the particle, plus energy of the emitted radiation) to be conserved. According 
to this view, the models so far proposed (GRW, CSL, QMUPL) are first approximations of more 
realistic models of spontaneous wave function collapse, yet to be formulated. 



B. The spectrum of the emitted radiation 

By using Eq. (|35p . with G^(k,t) given by Eq. (|4ip . we obtain the following expression for the 
time derivative of the emitted spectrum: 



d . . Xhe 2 1 

[<->col(,k, /X, tj 



dt colv ' r "' ' 167r 3 eo uJk 

2 



2m 2 + f3 2 Lu( 2(5 
— T^Tt — 2 i Ql 2\ ^ 27 — 2 i Ql 2\ L ° k sm 

uJk cos Wfct . (46) 



Since all observations are made over a period of time [10] much longer than the characteristic 
photon's frequencies, the two oscillating terms in the above expression average to 0. We are then 
left with only the first expression within brackets. 

The physically interesting quantity is the spontaneous photon-emission rate dT^/dk per unit 
photon momentum. This is obtained from dS col (\t, (J,, t)/dt by summing over the polarization states 
and integrating over all directions in the photon's momentum space. The final result is: 

d Tk _ Xhe 2 2 + {(5ck/m) 2 ^ 

dk 27r 2 eom 2 c 3 k 1 + {(5ck/m) 2 

It reassembles Eq. (21) of Q] (and Eq. (3.14) of Q|), when replacing eo — > 1/4-7T because of 
the different system of units used, and when taking A = (m/m N ) 2 Xo (m N is the nucleon mass) 



as assumed in the mass-dependent CSL model 



a. 



The only difference is the extra factor [2 + 



IS 



{[3ck/m) 2 \/[l + (f3ck/m) 2 ], the (3 dependence in which comes about because the result of 
has been carried out only to first perturbative order, while our result is exact (within the dipole 
approximation). For an electron, (f3ck/m) 2 ~ (9.47 x 10 _6 £'fc/KeV) 2 , where = hck is the 



energy of a photon of momentum k. Table I of [8| reports data from photons in an energy range 
between 11 and 501 KeV: our calculation shows that, in this range, the first-order perturbation 
theory is extremely accurate. 

Since Eq. (|47p is valid for finite times, it provides a trustable understanding of the radiation 
process within the limits of the dipole approximation, i.e. as long as the particle does not move 
too fast, or as long as the photon's momentum is not too large. By keeping only the leading 
terms in the relevant parameters, i.e. by setting j3 = 0, Eq. ([17]) reduces to twice the large-time, 
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first-order CSL expression of [l0(] and [llj. However, according to the argument of Sec. [H the CSL 
and QMUPL models should agree for sufficiently well localized systems 6 (with respect to the scale 
set by rc — 10~ 5 cm); the origin of this discrepancy will be the subject of further exploration. 

As a last comment, we note that Eq. (|47p predicts an infinite amount of energy to be emitted 
per unit time, as dT k /dk is of order 1/k for large k. This ultraviolet catastrophe is a consequence 
of the dipole approximation. One of the effects of the term e ik x in Eq. ([9]) is to temper the 
electromagnetic coupling for high frequencies; by replacing e lkx with 1, this effect is neglected. 
Accordingly, Eq. (|47p is not trustable anymore in the very large k limit. 



IV. THE HARMONIC OSCILLATOR 



When the particle is bounded by a linear force, the emitted spectrum takes a quite different 
expression. By inserting Eqs. (133 j) . and ignoring the term i = 1 which gives a runaway solution, 
one finds: 

ft 



I 

Jo 



Gi(k,t-s)Gj(k,t-s) 



£ 



(z£ + z t ,)(zf> - iu k )(z£> + iuJk) 



i=2 
3 



U k Z£ 



e=2 



H(-ito k )(z£ - iu k ) 2 

UkZe 

H(+iu k )(z£ + iu k ) 2 



Z — Z£ 
Z — Z£ 

~H{z) 



Z — Z£ 



z=z e 



Z — Z(i 

~H{z) 



(z e +z t ,)t 



z=z e 



Z = Z( 



D (z e +iu} k )t 



+ - 



LOT 



t. 



(48) 



H(-iuj k )H(+iui k ) 

The formula is rather cumbersome. However the terms in the first three lines contain exponentially 
decaying terms, which vanish very rapidly with time. For example — with reference to Eq. (|3L|) — 
the decay time is about 2.93 x 10~ 47 s for a 11 Kev photon. Accordingly, in the large time limit 
we have for the differential photon emission rate dT k /dk (where, as in the free particle case, we 
have differentiated Eq. (|35p over time, summed over the polarization states, and integrated over 
all directions in photon's momentum space) the following simple large-time expression: 

dV k A/ice 2 A; 3 



dk 



2vr 2 e m 2 {ujl - c 2 k 2 ) 2 + /5 2 c 6 /fc 6 ' 



(49) 



6 One can argue that the free particle case contradicts this assumption, as the wave function of a free particle 
rapidly spreads out in space; however, at least for sufficiently short times the approximation is correct. 
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Two comments are at order. The first important thing one notices is that Eq. (|49|) does not reduce 
to (|47p in the free particle limit. The reason for this incongruence can be traced back to Eq. (I48p . 
according to which the free particle limit (wo — ► 0) and the large time limit (t — > +00) do not 
commute, as one can prove by direct calculation. Prom the physical point of view, the reason for 
the discrepancy is that, in the large time limit, the particle has the chance to move far enough to 
feel the edges of the harmonic potential, no matter how weak the potential is. This mean that the 
particle is never really free, even in the limit ojq — > 0. As a further proof of this statement, one can 
note that by taking the free-particle limit at finite times, one does indeed recover Eq. (|47p . 

As a second observation, one can see that in the lowest order in the relevant parameters (/? = 0), 
the emission rate given by Eq. (|49p is of order 1/k for ck 3> ujq. This is reminiscent of the free 
particle case. However, the exact expression is of order 1/A; 3 , and the total emission is finite, 
contrary to what implied by the free-particle expression. The physical reason is that the binding 
potential works against the emission of high-energy photons, as the term e* k x in Eq. ([9]), which is 
neglected by the dipole approximation, does. 

The third relevant observation is that Eq. (|49l) shows a resonant behavior in correspondence to 
the natural frequency ujq of the oscillator. Indeed the peak of the resonance is very high, due to 
the fact that (3 2 c 6 k e is a very small quantity (confront the small value of (3 given in Eq. (|30p ) for 
k = ujq/c, where ujq is a standard frequency such as that associated to the hydrogen atom. Indeed 
such a great resonance is incompatible with experimental data and, as such, it would disprove this 
model, for any significant value of the collapse parameter A. However, the large value of the peak is 
an artificial feature of the model. It emerges as a combination both of the fact that the the energy 
levels of the harmonic oscillator are equally spaced, and from the dipole approximation, according 
to which transitions are allowed only between to consecutive levels. In order words, what happens 
here is that the noise excites the particle to an higher energy level state; in the de-excitation 
process only photons of energy fvjjQ can be emitted. In a more realistic model, also photons with 
any energy uItwq should be emitted, and the spectrum would have a more articulated resonance 
structure, where the peaks are less pronounced. An accurate spectrum would then display several 
resonances. 

To conclude, our analysis shows that, in presence of a discrete spectrum (e.g. the hydrogen 
atom), the differential photon emission rate due to the collapse process should show a typical 
resonant behavior, which has not been depicted by previous analysis. Although it is reasonable to 
expect that these resonances are highly suppressed, it is worthwhile analyzing such a behavior for 
the CSL model, by generalizing the previous results of 
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llj to the low-frequency part of the 
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spectrum. 

V. CONCLUSIONS 

We have analyzed the electromagnetic properties of both a free particle and of a particle bounded 
by an harmonic potential, within the framework of collapse models. By choosing a particularly 
simple, yet physically meaningful, model of spontaneous wave function collapse, and under only 
the dipole approximation, we have been able to solve the equations of motion exactly. 

In the free particle case, we have found a counterintuitive result: the particle's kinetic energy 
steadily increases in time, and at the same time it spontaneously emits radiation at a constant 
rate. Although this is in principle possible, as long as no conflict with experimental data emerges, 
such a behavior suggests that collapse models should be modified in order the temper (or eliminate 
entirely) the evident violation of the energy conservation principle. 

We have also found some discrepancies between our formula and those previously derived, 
through a perturbative analysis. The origin of these differences is not clear yet, and will be further 
studied in the future. 

In the case of a particle confined by an harmonic potential, the spectrum is modified and a 
peak emerges, in correspondence to the natural frequency of the oscillator. This feature suggests 
that also in more realistic situations (e.g. atomic systems) the spectrum should have a resonant 
structure, which is worthwhile analyzing. 

These results show that further analysis is required in order to better understand the electro- 
magnetic properties of charged particles in the CSL model. This is important both for clarifying 
the theoretical picture offered by collapse models, and also in the light of future experimental tests. 
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APPENDIX A: DERIVATION OF EQ. (31) 

The zeros of H(z) defined by Eq. ([301) correspond to the solution of the cubic equation: 

K m 

y ai = ' a2 = "^ 

By defining 



z° + a 2 z +aiz + a = 0, a = — — , ai = 0, a 2 = — — . (Al) 



q= \ ai ~\°^ r = ^ ( a i fl 2 - 3a o) - ^a|, si, 2 = \fr±y/q 3 + r 2 , (A2) 

the three roots can be written as: 

zi = (si + s 2 ) - ^a 2 , ^2,3 = (si + s 2 ) - ^ 02 ± i-^ (si - s 2 ) . (A3) 

This is the standard Cardan's method for finding the roots. In our case: 

3 2 1 2 A 4 4\ 1 k 2 A 4 m 2 \ 1 wgm 4 



if ujq <C 2m/\/27p\ as we have originally assumed. Then ^/g 3 + r 2 ~ uiqui 2 / \/27 (3 2 . Working under 
the same approximation we have: 

r ± yV + r 2 ~ __ M ± V27 — + — A5 

27 /3 3 V m 2 m 2 J K 1 



and 



v 3 p V ml m z 



From the above expression and from Eq. (|A3p . the approximate values of the roots given in ([31 
can be immediately derived. 
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